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Abstract 

In this work the induced van der Waals interaction between a pair 
of neutral atoms or molecules is considered by use of a statistical me- 
chanical method. Commonly this interaction is obtained by standard 
quantum mechanical perturbation theory to second order. However, 
the latter is restricted to electrostatic interactions between charges and 
dipole moments. So with radiating dipole-dipole interaction where re- 
tardation effects are important for large separations of the particles, 
other methods are needed, and the resulting induced interaction is 
the Casimir-Polder interaction usually obtained by field theory. It 
can also be evaluated, however, by a statistical mechanical method 
that utilizes the path integral representation. We here show explic- 
itly by use of the statistical mechanical method the equivalence of the 
Casimir-Polder and van der Waals interactions to leading order for 
short separations where retardation effects can be neglected. Physi- 
cally this is well known, but in our opinion the mathematics of this 
transition process is not so obvious. The evaluations needed mean 
a transform of the statistical mechanical free energy expression to a 
form that can be identified with second order perturbation theory. In 
recent works [Il0ye 2010] the Casimir-Polder or Casimir energy has 
been added as a correction to calculations of systems like the electron 
clouds of molecules. The equivalence to van der Waals interactions 
to leading order indicates that the added Casimir energy will improve 
the accuracy of calculated molecular energies. We here also give nu- 
merical estimates of this energy including analysis and estimates for 
the uniform electron gas. 
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1 Introduction 



The well known Casimir effect is commonly regarded to be due to ffuctuations 
of the quantized electromagnetic field in vacuum and has been widely studied 
in this respect. However, H0ye and Brevik considered this in a different 
way by regarding the problem as a statistical mechanical one of interacting 
fluctuating dipole moments of polarizable particles. In this way the Casimir- 
Polder force pLj between a pair of polarizable point particles was recovered 
P]. To do so the path integral formulation of quantized particle systems 
was utilized j3]. With the path integral the quantum mechanical problem 
at thermal equilibrium becomes equivalent to a classical polymer problem 
in four dimensions where imaginary time is the fourth dimension. The role 
of the electromagnetic field is then to mediate the pair interaction between 
polarizable particles. Later this type of evaluation was generalized to a pair of 
parallell plates, and the well known Lifshitz result was recovered [1] . Similar 
evaluations were performed for other situations [3 E]. Earlier this method 
was fruitfully utilized for a polarizable fluid [7]. 

The statistical mechanical approach opened new perspectives for evalua- 
tions of the Casimir force. Instead of focusing upon the quantization of the 
electromagnetic field itself one can regard the problem as one of polarizable 
particles interacting via the electromagnetic field. It was found that these 
two viewpoints are equivalent [21 IH El E] ■ 

The source of the Casimir free energy is the radiating electromagnetic 
interaction. For large separations of a pair of atoms or molecules retardation 
effects are important. Thus methods of field theory have been widely used to 
obtain the resulting Casimir force or energy for various situations [H]. Such 
situations are a pair of atoms, an atom and a wall, or two walls separated 
by some distance. Further the Casimir energy influences the shift of energy 
levels of atoms, and it is closely related to the Lamb shift as they both are 
consequences of the vacuum fluctuations of the electromagnetic field [9]. 

For short separations, however, retardation effects can be neglected, and 
on physical grounds it is well known that the Casimir interaction becomes the 
van der Waals interaction in this case. The usual way to compute the latter 
is to apply second order perturbation theory based upon the Schrodinger 
equation. 

Although the equivalence or the transition from Casimir to van der Waals 
interaction is well known from the physical point of view, it is not obvious 
or trivial from the mathematical side as the methods used are very different. 
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To obtain the Casimir interaction summation over imaginary Matsubara fre- 
quences are performed and the frequency dependent polarizabihty is needed. 
On the other side the van der Waals interaction can be expressed directly 
in terms of the matrix elements and eigenvalues of second order perturba- 
tion theory. Thus in this work we will show explicitly the equivalence of 
the Casimir and van der Waals interactions in the case of static interactions. 
To do so the statistical mechanical free energy expression for the Casimir 
energy to leading order is transformed such that it can be identified with 
second order perturbation theory. The equivalence obtained coincides with 
the interpretation that the Casimir force can be related to fluctuating dipole 
moments. The statistical mechanical method introduced in Ref. 0| is based 
upon this interpretation realizing that the electromagnetic field can be fully 
replaced by pair interactions between dipolar moments [H |5l El [8] . 

The mathematical equivalence of the Casimir and van der Waals interac- 
tions, which we will show below, gives a close and direct connection between 
the Casimir theory and standard quantum mechanical perturbation theory. 
This means that the Casimir interaction can be interpreted as a perturbing 
correction upon quantized many-body systems of interacting particles. Thus 
it is of interest to take this energy into account for a more general situation 
of such systems like electrons in molecules. This direct and close connec- 
tion seems not to have been noticed previously. Radiation effects, if taken 
into account, will be minor corrections. Recently work in this direction was 
initiated by Buenzli and Martin pDJ and by Jancovici and Samaj [llj. 

Earlier the path integral formulation of fermion and boson fluids was 
considered by Chandler and Wolynes [12] and by H0ye and Stell [13]. But 
the latter investigations were restricted to uniform density and instantaneous 
interactions. Present standard methods to evaluate many-electron systems 
utilize Hartree-Fock (HF) and density functional theory (DFT) [S] . By these 
latter methods exchange effects are taken into account, and approximations 
for local correlations are made. In separate works we have recently reconsid- 
ered this problem by adding the corresponding Casimir interaction to ab ini- 
tio calculations of interacting many-body systems like electrons of molecules 
[T^ [TB] . In Ref. [IE] we also were able to include the radiating electromag- 
netic field, and it was found that an ionic fluid of electrons bound to atomic 
nuclei can be identified with a dielectric one. 

With the Casimir energy included non-local correlations between elec- 
trons are accounted for. These correlations correspond to the induced ones 
between the fluctuating dipole moments (or polarizabihty) of the pair of point 
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particles considered in the this work. To estimate the influence of these corre- 
lations we will make crude estimates for atoms at close contact, and further 
make an analytic approximation with numerical estimates for the uniform 
electron gas. 

To compare the statistical mechanical method with standard second order 
perturbation theory some obvious results of the latter are needed and are 
therefore established below. Thus in Sec. [2] we write down the well known 
expression for the static polarizability due to a particle or electron bond by a 
potential. This is expressed in a standard way by means of energy eigenvalues 
and matrix elements between eigenstates. 

In Sec. [3] the appropriate pair correlation function in imaginary time 
needed for the statistical mechanical method is obtained. To do so the corre- 
lation function in imaginary time is established in view of the path integral, 
and Kubo relations are used. With the simplifled dipolar type interaction 
used in Sec. H] one flnds that its Fourier transform in imaginary time is the 
well-known dynamic polarizability. 

In Sec. mthe statistical mechanical method is used to obtain the Casimir 
interaction between a pair of particles. This is the result we want to com- 
pare with standard quantum mechanical perturbation theory to second order 
for static interactions. To do so the leading term in the Casimir interaction 
is considered and explicit summations with respect to the Matsubara fre- 
quencies are performed. Then the expressions obtained for the dynamical 
polarizability are used, and we show that the Casimir energy can be ex- 
pressed in terms of the usual matrix elements of quantum mechanics. The 
result depends upon temperature in general. 

In Sec. [5] the statistical mechanical result obtained in Sec. IHis verifled to 
be equivalent to standard quantum mechanical perturbation theory to second 
order. 

In Sec. [6] the situation is generalized in a straightforward way to more 
particles or electrons bond to a potential. Then properties of fermions (or 
bosons) have to be taken into account to obtain the appropriate correlation 
function in imaginary time. Again this is found to be the corresponding 
dynamical polarizability. 

In Sec. [7] the statistical mechanical approach is again identifled with sec- 
ond order perturbation theory for the more general situation with several 
electrons, which is most relevant for realistic many-body systems. 

In Sec. [8] numerical estimates of the Casimir energy for neighboring pairs 
of atoms in molecules due to induced correlations are made. Then a pair 
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of polarizable particles at close contact, but otherwise regarded as point 
particles, are considered. 

In Sec. [9] the homogeneous electron gas is considered and an analytic 
estimate of the correction to ab initio calculations are made. An explicit 
approximation that covers all electron densities is obtained. Especially for 
low density it is found that the ground state energy of plasma oscillations 
is the main contribution besides the binding energy. One may speculate 
whether this has relation to high-Tc super-conductivity that is considered to 
be connected to properties of electrons at low density. 

In Sec. [To] numerical estimates of the analytic result obtained in Sec. [9] 
are made. One estimate is for the density of conduction electrons in a metal. 
In this connection it is noted that the Casimir energy can be associated with 
the random phase approximation (RPA) for quantized many-body systems. 
Work has been performed on the RPA to improve it further in view of best 
available "exact" results. So the analytic approximation of Sec. [9] is com- 
pared with one such result for high electron density. Further the statistical 
mechanical viewpoint indicates new ways to improve the RPA. This will be 
discussed shortly. 



2 Static polarizability 

The statistical mechanical method requires the appropriate correlation func- 
tion of a reference system as input. In the present case with pair interactions 
of dipolar type this will be the dynamical polarizability as derived in the next 
section. First we will recover the simple result for the static polarizability. 
So consider a particle in an external potential like an electron bound to an 
atomic nucleus. It is perturbed by an electric field along the x-axis. This 
contributes 

H' = -qx£^ (1) 

to the Hamiltonian where q is charge and x is position. The applied field 
Sx produces a net dipole moment Px in the same direction provided the 
unperturbed system possesses spherical symmetry. By standard perturbation 
theory one easily obtains the static polarizability a either via the average 
position (x) as Px = q{x) = aSx or via the perturbed energy to second order 
given by E2 = —\Px£x- In either case one finds [IT] 

« = 2Ef^ (2) 



mO 
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provided the particle is located in the ground state as is the case for low 
temperatures, T — )■ 0. Here we have introduced the abbreviations 

Mynn = {4>m\x\(pn) {(Pn\x\4>m) 

where the 0^ = <Pm{x) are the unperturbed wavefunctions and are the 
corresponding eigenvalues for the energies. The 0o and Eq are ground state 
quantities. It can be noted that in this work all states are considered to 
be bound states of the potential to avoid possible problems connected to 
unbound states with continuous spectrum. 



3 Dynamic polarizability 

With the path integral representation the correlation function of interest can 
be obtained from the one particle correlation function in imaginary time A. 
This function is given by [151 fT6] 

G{X,Xi,X2) = ^Y.^rnn{>^,Xl,X2) (4) 
mn 

where 

Fmn{\, xi, X2) = 0™(xi)e-(''-^)^™0;;(x2)0„(x2)e-^^"0:(xi) 

and Z =^e-^^'= (5) 

k 

recognizing that in the path integral picture the quantized particle represents 
a "classical" polymer in the imaginary time direction. In this direction it has 
an extension /3 = 1/(A;bT) where T is temperature and ks is Boltzmann's 
constant. The relation between A and real time tist= —ihX. (For simplicity 
one dimensional notation is used here and below. By a mistake the A was 
defined with opposite sign in Ref. 0, but this did not influence previous 
applications where (apart from a constant) the operators A and B of Eq. ([9]) 
below were equal. The mistake was corrected in Appendix of Ref. jl6j.) 

The dipole moment and thus a is now found from the appropriate 
response function via the Kubo relations |2l |T6l |18]. Then we need the 
average 

,(A) = g^x,X2) = g^Jx,X2GiX,x„X2)dx,dx2 

= 5]p„M„„e^^- (6) 
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where 

Pm = e-^^-/Z (Ep™ = 1)- (7) 



m 



The Prn is the probabihty for the particle to occupy energy level m at thermal 
equlibrium. Further the Fourier transform is needed 

? 1 

~g{K) = g{\)e'^^d\ = Y,{pr.-Pm)M^^——— 

^ mn ^ ^mn + f'J^ 

where the symmetry properties M^n = and fi^n = ~^nm have been 

used. The K = 27rn/P with n integer are the Matsubara frequencies. 

The resulting dipole moment Px = q{x) is now related to the applied 
(time varying) electric field £x via the response function [21 [161 H] 



^B^(t) = ^Tr(p[A,i?(t)]) (9) 

where p is the canonical equilibrium density matrix p = e~^^ /Z {Z = 
Tr{e~^^)), and the B{t) in the commutator is the time-dependent operator 
B[t) = e^H/hj^^-\H/h ^i^Qj-e X = it/Ti. For a time- dependent perturbation 
H = —AF{t) the change in the quantity B is then given by 

t 

{AB{t))= J ^BA{t-t)F{t)dt. (10) 

— oo 

As noted above the A by a mistake was defined with opposite sign in Ref. [2], 
but it was corrected in Ref. [16]. By this correction the Schrodinger equa- 
tion in imaginary time transforms to an equation that describes diffusion 
for increasing A. With this the correlation function in imaginary time and 
response function in real time are 

g{\) = Tr(pB(t)A) 

VBAit) = ^W + X)-g{X)] (11) 

As found in Appendix B in Ref. [2] the Fourier transform of the correla- 
tion function in imaginary time equals the Fourier transform of the response 
function in real time, i.e. 

ifBAicu) = ~g{K) (12) 
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in the common region where lm{u)) < and |Im(ii')| < C (C > 0) with 
K = ihu. 

In the present case with perturbation ([1]) and polarization = q{x) = 
(AB) we have the operators A = qx = qx{0) and B(t) = qx{t) by which the 
correlation function will be the one given by Eq. ([6]) with Fourier transform 
OH]). With Pxi^jj) = Oi{K)£,j.{(jj) {K = ihu) it thus follows that the frequency 
dependent polarizability is given by 

a = a{K) = g{K). (13) 

For K = and zero temperature the static result is recovered. 

4 Induced interaction between a pair of par- 
ticles 

Consider a pair of particles located at different positions in space each with 
polarizability given by f[T^ . Then assume there is a mutual perturbing in- 
teraction between these particles given by 

$ = ifjq^XaXb (14) 

where ijj = ijj{r) is a quantity that depends upon spatial separation while 
Xa and Xh are internal coordinates of the two separate binding potentials. 
Interaction (fT4|) is similar to the dipole-dipole interaction where qxa and 
qxb correspond to dipole moments restricted to the x-direction. We will use 
the one dimensional interaction ( !T^ instead of the full dipolar interaction 
to simplify the derivations below a bit. The purpose is to transform the 
statistical mechanical free energy expression to a form that can be identified 
with second order perturbation theory. 

The free energy change AF due to the interaction can be found from |2] 

= ^ E 1^(1 - M)') ^ E(«^)' (15) 

to second order in a = a{K). For K = expression f|T5l) is just the classical 
expression for the free energy change due to a pair of classical harmonic oscil- 
lators that interact. This is Eq. (3.4) of Ref. [2]. In the quantum mechanical 
case one has to sum over the Matsubara frequencies K = 27Tn/P {n is inte- 
ger) as follows from the path integral. More generally the interaction can be 
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time-dependent like the radiating dipole interaction. Then the Fourier trans- 
form ip{K) with K = ihu is needed. However, to enable comparison with 
standard quantum mechanical perturbation theory the statistical mechanical 
expression should be restricted to the static case, i.e. ip = ip{0) for all K. 

It can be noted that the term to second order in a in expression f lTSj) is 
consistent with the original derivation of the Casimir-Polder force [T] when 
modified such that the ip is replaced with the radiating dipolar interaction 
that depends upon orientation of dipolar moments and frequency too, as 
done in Ref. [2]. Casimir and Polder used the methods of quantum electro- 
dynamics to obtain this force. Wennerstrom et al. made a closer study of 
the temperature dependence of the Casimir interaction where both a and ijj 
depend upon K It can be noted that their expression (5) (except for 
the n = term that is taken out) for this interaction between a pair of par- 
ticles is the same as expression (5.15) obtained by the statistical mechanical 
method in Ref. [2]. (A factor 2 is missing in the denominator of Eq. (5.15) of 
the reference, but is present in its Eq. (5.14).) This is seen when the r and 
ax of Ref. [2] are identified with the nxn and a(0)/(l + {An)"^) of Ref. [12] 
where A = 2'KkBT / {Tiujq).) 

In the present work we want to verify that the Casimir interaction for 
the specific case of static interactions, i. e. when ip does not depend upon 

is identical to the van der Waals interaction that follows from standard 
quantum mechanical perturbation theory to second order. Then expression 
f ll5p must be transformed to a form that is recognizable and can be identified 
with the latter. This is done by explicit summations in f lTSj) with respect to 
the Matsubara frequencies K = 2Tm/ (5. We need the sum {K = 2nn//3) 

I3^Q? + K^ 2sinh(/3fi/2)' ^ ' 

It follows from the backward transform of expression ([8]) for A = noting 
that the result should be the average of (5^(0) + g{f3))/2 since implicitly the 
g{X) is assumed periodic. With this we obtain 



1 ^1 ^2 ^i^l2 1 



1 



1 



- ni] 



cosh(/3^]l/2) cosh(/3f)2/2) 
^^^^^Vk/2)~ %inh(/3f]2/2) 



4sinh(/3fii/2) sinh(/3^]2/2) 



;i7) 
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X 



sinh(/5((]i + n2)/2) sinh(/3((^i - VL2) /2) 



In the limit /3 — )■ 00 this result simplifies to l/[2(|f2i| + |f^2|)]- 

By use of expressions ([8]) and (fT7|) we now find for the free energy contri- 
bution (ITSl) 

= -^Eiam' = -"i^Y.T.^rnuMkliPmPnPmY" 



K 



ran kl 



smh{(5{yLmn + ^ki)_/])_ ^ sinh(/3(a 



kl 



mn kl 



PnPl - PmPk _^ PnPk - PmPf 



^mn ^ 



-i^'^miPnPi 



M^nMkl 



mn kl ^mn ^kl 



(18) 



Here Pn—Pm = '^{PnPmY^'^ sinh(/3i7mn/2) has been used {pn = e~^^"/Z). The 
last equality in Eq. ( |T8l) follows by use of the symmetries Mmn = Mnm and 

^mn = — ^nm = Em — En hj SUmmatioUS. 



5 Second order perturbation theory 

For comparison the expressions of standard quantum mechanical perturba- 
tion theory are needed. For two separate particles with interaction f |T4l) the 
perturbed energy of the eigenstates n and / is 

Ai?nz = -il'^r E , = -r E 7^^^^ (19) 

mk -^"i + Ek — [En + El) \lnm + "fc« 

to second order. Here {mk\xaXb\nl) = {m\xa\n){k\xb\l) is used. (The contri- 
bution to first order with m = n and k = I i assumed zero.) For T = it is 
easily seen that the free energy flTK]) coincides with expression flTI?]) for Aii^oo 
as then the system is in its ground state with n = / = 0. 

For finite /3 one should consider the resulting partition function 



nl nl 

(^1 -^Ee"'^''"^'''^^^^^ + •••)• 
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With free energy F = — In Z^j /3 the perturbed part of it will be 

^^ = ^Ee-''^''"'''''^A^n^ + ---- (21) 

nl 



where expression f[T9|) is to be inserted. One sees that the statistical mechan- 
ical result f lTS]) coincides with this equation. This verifies the mathematical 
equivalence of the two methods for static interactions. 

6 Polarizability with several particles 

With several particles in a common potential, e.g. electrons in an atom, one 
has to take into account the symmetry properties of fermions (or bosons). 
However, the eigenstates can again be represented (approximately at least) 
by one-particle wave functions. But the correlation function (jl]) is modified 
into 

G'(A,a;i,a;2) = CY^ jnfniFnin{\xx,X'2). (22) 

mn 

Here 

= 1 + l,-,E. C^=e^^ (23) 

where /i is the chemical potential. (For bosons /„ = 1/(1 — ('e"^'^").) The 
/„ factors include all sums of path integral polymers that are tied together 
into coils. 

With fl22p it immediately follows that the expression for the polarizability 
as given by expression ([HD can be kept if p„ for fermions is redefined to be 

p.. = ACe-*'^" = (24) 

which is the probability for level n to be occupied by a particle. (Thus 
definition ( 12^ for p„ replaces the one by Eq. (J?!).) With this we again have 



a{K) = g{K) = E(Pn " ^"f • (25) 

mn ^ 'mn ^ 



noting that CUfn{e ^^"^ - e ^^"') = Pn - Pn 

It may be noted that expression fl25|) can be interpreted as a sum of 
independent contributions from each of the particles (or electrons) present. 
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This reflects the one-particle nature of the wave functions used. This is seen 
by again using Mmn = M^m and Q^n = -^nm by which the -p^ in Eq. <^ 
can be replaced with p„ to obtain 

(^{K)=2\^Pn ■ (26) 

mn mn ~r -i v 

Further the average number of particles that contribute to a{K) is N = 
J2nPn- Thus expression ( l26l) shows that all occupied levels can be considered 
to give independent contributions to the polarizability. 

Notably all terms in the sum (125|) are positive, and the most significant 
contributions come from energies close to the fermi energy (or chemical po- 
tential) fj, where Em > and ^ > En or vice versa such that fmPn ~ 1 and 
fimn is small. On the other hand the rearranged expression fl26p apparently 
obtains significant contributions from all levels below the fermi level. But 
since flmn can have both signs these contributions thus mostly must cancel. 



7 Induced interactions and perturbation with 
several particles 

With polarizability given by expression ( 125|) the results of Sec. |4]are by tak- 
ing the change of Pn into account, generalized in a straightforward way. With 
Eq. ([MD one now finds p„-pm = Pnfm-Pmfn = '^{fnPnfmPmf'^ sluh (/3fi^„/2) 
since Pn = (C/nJ5„)^/2e-^^"/2 and /„ = {(-^fnPnY^'^e^^-/^ Thus in the 
{PmPnPkPiY^'^ term of Eq. f|T8l) one can replace the pi {i = m,n,k,l) with 
fiPi, and with this result f|T8l) modifies into 

AF = -V^^ E E fmPnfm^^^^^. (27) 

Concerning second order perturbation theory the situation may be less obvi- 
ous. The problem is that one-particle eigenstates can be occupied by other 
particles. Thus for fermions one may assume that the state of interest is 
perturbed only via vacant levels. With this modification Eq. f|T9|) for the 
perturbed energy levels is changed into 

AE„z = -^^E/r^^#^- (28) 

mk ''mn ~r ' 'fci 
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For fermions the grand partition function for the particles in each of the 
potentials (or atoms), that here are assumed equal, will be 

Z = Y[{1 + Ce-^^") (29) 

n 

Thus the total grand partition function with perturbing interaction will be 

Z, = ]^(i + ^e-'^i^" + Ce-''^' + CV('^^"+^'+^^"')) 

■n,l 

= Z\l- PY,PnPl^Enl + ---) (30) 

n,l 

with pn given by (!24|) . With free energy —fi^^lnZt this again leads to the 
free energy contribution 

= Y,PnPi^E^i (31) 

n,l 

where now the p„ is given by Eq. f l2^ and AEni is given by Eq. fl28l) . With 
this one sees that the statistical mechanical result fl27|) is the same as result 



8 Numerical estimate for polarizable atoms 

As mentioned earlier the Casimir type energy gives a leading energy con- 
tribution to ab initio computations of of molecular energies. This type of 
contribution, which is a correction to mean field theory, was introduced for 
classical many-body systems with a perturbing interaction that was assumed 
weak and of long range where the inverse range 7 is the perturbing parame- 
ter [20j . For interacting harmonic oscillators it turns out that this correction 
alone can give the exact answer, so for this case expression (fT5|) (before ex- 
pansion in ip) is exact. 

Clearly, accurate results for the energy corrections to molecular energies 
will require extensive numerical evaluations on inhomogeneous systems. Thus 
we here can only give crude estimates. One way is to make an estimate of the 
van der Waals interaction between a pair of atoms at contact. With static 
dipolar interaction the T = result for a pair of harmonic oscillators is after 
integration with respect to K given by Eq. (1) of Ref. [T^] as 

V{R) = --huj,'^. (32) 
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The a(0) can be given a crude estimate from the susceptibihty of a polarizable 
hquid. For low dielectric constant e it is approximately given by (in Gaussian 
units) 

e-l^ 47rpa(0) = Ai^pa'^'^ (33) 

where a is the (hard core) diameter of an atom. An order of magnitude 
estimate for the susceptibihty may be to put e — 1 1.0 at (dimensionless) 
fluid density pa^ ~ 0.4 by which a(0)/(T^ ~ 0.20. Further one might assume 
the resonance frequency to give an energy quantum Tluq ~ 1.5 eV. With this 
the van der Waals energy (|32|) at contact R = a would be 

V(a) ^ -45meV. (34) 

A similar estimate can be obtained for instance on basis of the Lennard- Jones 
(LJ) interaction 0(r) = AeLj[{a / Rf^ -{a / Rf] of Ar (argon). The mmimum 
value of is (pmin = —^Lj- For Ar the critical temperature is Tc = 151 K, 
and for the LJ fluid tLj/lksTc) ^ 0.8. Thus with ks = 1.38 ■ 10'^^ 3 /K = 
8.63 ■ 10~^eV the attractive van der Waals part of the LJ interaction at 
contact R = a becomes 

- 4eLj ^ -42 meV. (35) 



9 Analytic estimate for an electron gas 

It will be more satisfying to base the estimate more directly upon the correc- 
tion to ab initio computations. But for the inhomogeneous fluid of electrons 
in molecules this does not seem feasible in a simple way. However, we find 
that an estimate can be made for the homogeneous fluid of electrons (in a 
neutralizing background). In view of 7-ordering in terms of the parameter 7 
mentioned above, the electrons without interactions form the reference sys- 
tem, which then is an ideal fermi gas. Since the electrons in the reference 
system or fluid are free to move the corresponding susceptibility e — 1 will 
be infinite at zero frequency (as for a metal) in this case; so there is reason 
to expect a larger correction. 

The correction AF to the free energy per unit volume due to the perturb- 
ing interaction is given by Eq. (5.8) of Ref. [12]. Now taking spin degeneracy 
g = 2 into account (both up and down spins contribute) this becomes 

- /3AF = / = X: / ln[l - gS{K, k){-m)]- (36) 
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The Fourier transform of the electrostatic Coulomb interaction in SI units is 

m = (37) 



eok^ 

where q is the unit charge of electrons and Eq is the permittivity of vacuum. 
In the uniform case the S{K, k) is the structure factor for the ideal gas of 
free fermions for spins either up or down. (Different spins are uncorrelated 
in the reference system.) It is given by Eq. (3.4) of Ref. ^Lbj or its modified 
version Eq. (63) of Ref. [16] which is 



S{K, k) 



c 



(2vr)= 



A 



X -Y 



with (m is particle mass) 



7^2 + A2 {l±CX){l±CY) 

{hkf 



(38) 



A = E{k") ~ E{k'), E{k) 



2m 



k" = k - k'. 



Y 



-l3E{k") 



(39) 



c 



where fi is the chemical potential (Fermi energy at T = 0). 

One may expand expression ( l36l) to separate out its leading term which 
is the exchange energy 

9l3 



S{0,k){-^lj{k))dk. 



2(27r)3 

As given by Eq. (3.1) in Ref. [I5] (for A = or A = f3) 



SiO,k) = -J2SiK,k) 



K 



X 



{27ifJ (i + cx)(i + cr) 



dk! 



(40) 



(41) 



which is the equal time correlation function. In the limit /3 — i- oo one finds 



i2n) 



dk' 



k'<kj 
k">kf 



i2n) 



k' 



dk' 



k'<kf 



k'<k 



f 



1 u -\ u 



4 
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gSiO,k) 



p for 2 < u with u = — 

kf 



k"<kj, 

for < -u < 2, 
k 



(42) 
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where kf is the value of k at the Fermi energy, and p is the total particle 
density 

p = -l- [ dk=4^kl (43) 
^ (27r)^ J 3(27r)3 ^ ^ ' 

k<kf 

The p-term in Eq. ( H2l) is a self-energy term ((5(r)-function in r-space) and 
should be deleted from Eq. fHOj) at least with the Coulomb interaction where 
it will give an infinite contribution. The remaining part of Eq. fH2|) is the 
known equal time pair correlation function of fermions at T = 0. (It is the 
same as the overlap volume of two spheres each of radius 1 with centers 
separated by a distance u.) Inserting this and expression (ET]) for ipi^k) gives 
the known exchange energy per particle {g = 2) 

_AFe. _ 47r 3q% _ ^{huj,f 
^'"^ p (27r)3 8£o 32p ^ > 

where Up is the plasma frequency and is the Fermi energy 

' meo 2m 



The higher order terms of Eq. (1361) give the correction to the Hartree-Fock 
result. In the classical case they give the well known debye-Hiickel theory 
of ionic solutions and plasmas where Debye shielding is important. For this 
situation gS{0, k) = (3gS{0, k) = f3p (since then (3, p, (X — )• 0). However, for 
T = the correction term turns out to change character. Then the classical 
term K = alone gives no contribution. Instead as noted in Ref. [TS] the 
resulting correlation function reflects the presence of plasma waves. Then 
the A in the denominator of Eq. ( |38|) was considered negligible and to a first 
approximation equal to zero. However, to avoid divergence the A can not 
be neglected here. But as a simplification the k will be regarded small, by 
which A can be approximated by 

A ^ . 46 

m 

To be able to handle Eq. (138|) in an analytic way we have to approximate 
further. Thus at T = {Ek -> iP/i^n)) J dK) 

1 r 2A 

' ^ ■dK=l, (47) 



2nJ + A2 
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and as simple approximation Eq. 

S{K, k) = 



may be written as 
A2 2S{0,k) 



K^ + A^ lAI 



(48) 



Used in Eq. (HTj) (for T = 0) this gives an identity for S{0, k) and is thus con- 
sistent independent of A. But for products of S{K, k) it can only serve as an 
approximation that may preserve the main features of the energy correction. 
So some approximation for A related to the integration with respect to k' is 
needed. We choose to take the average of A^ inside the fermi surface. Then 



/ 



((kk')^) = \e{k'^) = \ee 

from which with /i given by Eq. ( I45l) one gets 



A^ = {cijJLuY with ci 



0.8944, 



u 



k 



(49) 



(50) 



With ^(O, k) given by Eq. ( l42l) we with approximation f j48l) find (for small 

A^ 2p\u g2 



u) 

A{K) 
where 



-gS{K,k){-i;{k)) 



A' 



a 



K'^ + A2 ci/iu eofc^ K'^ + A2 u'^ 



hUr, 



a 



C2- 



C2 



'A 

4ci 



0.9157. 



(51) 



(52) 



Here it can be remarked that when plasma waves are considered, the A^ 
in the denominator can be neglected when compared to K^, and one finds 
A{K) = —Up/u"^ which would require C1C2 = 1 or ci = 4/3 in contrast to 
the value given by Eq. f l50|) . However, in our approximation we will keep the 
latter reasonable value for A . 

Integration of Eq. fl36|) with respect to K with the linear term taken out 
further gives 



2 



■1 + 1^ 



u 



u 



(53) 
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By final integration with respect to k the integral will diverge due to 
our approximation that is limited to small k. Thus the integration should 
be cut, and here the cut will be defined by regarding the low density limit 
where the ground state energy of plasma oscillations hup/2 will dominate 
besides the binding energy. The cut with respect to or u will be made at 
u = uq such that each particle gets this ground state energy. For a large 
the square root term of Eq. flS51) times A is then A{a/u) = ciC2hujp. With 
(1/(27?)^) (ik = {3/g)pu'^ du (with u = k/kf and g = 2) this implies 



3 7 2 , 1 s 
1 = - / C1C2U du = -C1C2M0, Uq = 1.3466. 

9-1 9 



(54) 



The correction to the free energy, where the linear term is separated out 
from Eq. fl36|) . now becomes 



fc 



AF - AFe, 
P 



Uo 



9 



fc{k)u^ du. 



(55) 



To integrate it is convenient to introduce a new variable of integration x 
u/a with upper limit of integration 



Uq 

Xo = — . 

a 



With this and relation dHU) one has 



(56) 



Ci/ia = Ci/xa 



CiC2xl 



Xq 



(57) 



so with fc{k) given by Eq. (153|) we obtain 



fc = -hur, - - 
JC 2^8 




Xq 



Xn 



2xo 

Xo 



■ (58) 



This result is verified by putting xq = x in expression (1581) and then differ- 
entiate x^fc/xo with respect to x to obtain {3/g)a^x'^ fc{k). 

It can be noted that result fl58|) is a valid approximation for any density. 
For high densities (small a) the energy of the electron gas is dominated by 
the kinetic energy of free fermions and the exchange energy. For low density, 
however, where a — )■ 00 [xo — )■ 0) this is no longer the case. Then the ground 
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state energy of plasma oscillations in addition to the binding energy becomes 
important. The latter is the exchange energy plus a similar contribution 
from the I/xq term at the end of Eq. fl58|) . that may be understood as a 
contribution from induced correlations. 

The above result may have relevance for high-Tc superconductivity where 
a low density strongly correlated electron gas is expected to play a role. 
In view of Eq. flSSl) strong correlations mean strong binding energy plus 
ground state energy of plasma oscillations while the Fermi energy becomes 
less significant. Due to the dominance of the binding energy the electrons 
may prefer to form a regular rigid lattice. This lattice of electrons may adapt 
to and thus stick to the underlying atomic lattice structure, or it may not 
do so, depending upon whether net energy is lowered or not. In the first 
situation the system may be an isolator while in the second case the electron 
lattice may be free to move relative to the atomic lattice. If so, the system 
becomes superconducting. However, these speculations, if relevant, will need 
further investigations. 



10 Numerical estimate for the electron gas 

As a numerical example we will consider the energy correction for the electron 
gas corresponding to one conduction electron per atom in copper (Cu). With 
atomic weight 63.54 , mass density 8.96 ■ lO^kg/m'^, and Avogadro's number 
Na = 6.022 ■ lO^'^mol"^ the particle density becomes 

p = 8.5 ■ lO^^m"^ (59) 

Further with electron mass m = 9.11 ■ 10"'^^ kg, electron charge q = 1.602 • 
10-^9 As, eo = 8.854 ■ lO'^^ (As/(Vm)), and h = 1.054 ■ lO'^^ Js one finds 
(leV= 1.602- 10-19 J) 

^Wp = 10.82 eV, /i = 7.04eV. (60) 
Thus with Eqs. ([56]), §M>, and one gets 

xo = = 0.957. (61) 

C2 nup 

When this is inserted in Eq. f l58|) the correction or exchange-correlation en- 
ergy becomes 

/,^-1.4eV. (62) 
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As a comparison the average particle energy of the reference system free 
electron gas is (3/5)/i = 4.22 eV while the exchange energy fl44l) is fex = 
— 4.68eV. Thus for the uniform electron gas the estimated magnitude of the 
correction beyond the exchange energy is quite significant for the density of 
conduction electrons in metals. In view of result f l58|) this significance will 
increase for lower densities of the electrons (at T = 0). 

It here can be noted that when comparing the results of Ref. [12] with 
other approaches one finds that the Casimir energy, at least when electro- 
static interactions are used, is equivalent to the results of the random phase 
approximation (RPA) for quantized many-body systems [2T| 122] . With the 
RPA it is expected that the energy functional should be optimized or made 
self-consistent with respect to its RPA part too. This results in a formidable 
computational task which limits the method to simple and small systems 
[22] [23] . However, from the viewpoint of Casimir energy used in this work 
there should be no need for such optimization of the RPA part since from the 
statistical mechanical derivation and 7-ordering |20] this part is the leading 
correction to HF and DFT computations. Its influence upon wave functions 
will be similar. At optimum a change in wave functions will to linear order 
in this change have no influence upon the energy; it will be quadratic and by 
that of higher order. Thus it can be neglected, and the formidable compu- 
tational task just mentioned can be avoided. Likewise no such optimization 
were needed to obtain the Casimir energy where possible non-linear parts 
of the dielectric constant are neglected. This is clearly the case for parallell 
plates filled with electrons where change in charge density at the surfaces 
due to electron interactions were neglected |23]. Thus in this respect similar 
changes in molecules can be neglected. 

The RPA is the quantized version of the Debye-Hiickel theory for ionic 
systems. As pointed out by Lein et al. the RPA approximation gives too low 
energy |23j. (The situation is similar for the classical Debye-Hiickel theory.) 
This they correct by modifications of the correlations at short range. Related 
modifications in view of the statistical mechanical picture are indicated at 
the end of Sec. 5 of Ref. [TJ]. Especially we here have in mind the simple 
properties of the direct correlation function that have been experienced for 
classical systems. These properties are expected to carry over to the "poly- 
mer" path integral of quantized systems as has been demonstrated for hard 
core classical polymers [25] . Then the hard core condition defines the direct 
correlation function or an effective interaction inside the hard core diameter 
via the Ornstein-Zernike integral equation. The perturbed correlation func- 
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tion of electrons (for the mixture of spins up and down) should also fulfill 
a hard core condition for small r (exact for r = 0). The /^.^(g,^) function 
(or the corresponding local field factor G{q,u))) of Ref. [23] serves a similar 
purpose to obtain an effective interaction, and various approximations to it 
have been considered. As far as we can see they do not include the use of 
a hard core or an effective hard core condition which also can be applied in 
the non-uniform case. 

In Ref. [23j| the various corrections to the RPA were compared with "ex- 
act" results. One such case is for = 1 with 

^{Tsaof = - (63) 
o p 

where = A-ne^h'^ / {mq^) = 5.29 ■ 10~^^m is the Bohr radius. This gives 
particle density 

p = 1.61 ■ 10^°m-^ (64) 

for which we find 

hujp = 47.1 eV, /i = 50.1 eV, xq = 1.564 (65) 

with correction or exchange-correlation energy 

/c^-2.4eV. (66) 

Various "exact" results referred to in Ref. [23] are -1.63 eV, -1.62 eV, -1.53 eV, 
and -1.64 eV. In view of the approximations made, result fl66|) is reasonable. 

However, the approximate simple expression fl53|) can be a basis for more 
accurate evaluations. In this respect we have made a comparison with Fig. 2 
of Ref. [23] where the energy distribution edq) = 3fc{k)u'^ (g = k) is plotted 
as function of q/{2kf) = u/2 for the case = 4. For small m — )■ our fc 
follows both the RPA and "exact" curves, while for larger u it becomes more 
negative as reflected in result (!66|) . Then to obtain the RPA curve one in 
approximation 0481) can replace the S{0, k) with an effective quantity (as a 
substitute for exact evaluations). Further to obtain the "exact" curve from 
the RPA curve the interaction fj37|) is modified for large k. This modification 
is a small r contribution to the direct correlation function (which indirectly 
defines the fxc{,k,(jj) mentioned above). So multiplying the fj37|) with the 
simple factor 1/(1 + ck'^Y (c =const.) already gives a good approximation 
to the exact curve. The challenge and numerical task will be to predict 
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quantities like the constant c on basis of the hard core condition mentioned 
above. The predictive power of the hard core condition will be especially 
useful for the non- uniform case of electrons forming molecules where "exact" 
results for comparisons are less available. 

11 Summary 

We have evaluated the induced van der Waals interaction between two point 
particles (e.g. two atoms) by a statistical mechanical method. Comparing 
with standard quantum mechanical perturbation theory to second order the 
results are the same. This seems obvious from the physical side, but math- 
ematically this is not so. The equivalence is restricted to short separations 
where retardation effects can be neglected. 

The statistical method can also be applied when retardation effects are 
present, and the resulting interaction is the Casimir-Polder one. The equiva- 
lence shows that the Casmir interaction as obtained by the latter method is 
an effect relevant for quantum mechanical evaluations of many-body systems 
like many-electron molecular systems. The leading term in the expansion of 
Eq. f|T5|) corresponds to second order perturbation theory. In addition there 
are higher order terms, and time dependent interactions like the electromag- 
netic one can be included. Thus in recent works we added the corresponding 
perturbation to ab initio calculations of systems like the electron clouds of 
molecules [151 HH] • This perturbation gives the leading energy correction to 
the well-known Hartree-Fock and density functional theory [H] . In this way 
induced energies due to non-local correlations are accounted for, and shield- 
ing of electric charges, which is standard Debye-shielding in the classical case, 
is embedded when contributions corresponding to all terms of expression (ITSj) 
are included. 

To see the magnitude of this perturbation we have estimated the inter- 
atomic induced energy for a pair of atoms at close contact. On basis of 
their van der Waals interaction this energy was found to be of the order of 
—40 meV. Further we have analyzed the equations for the uniform electron 
gas and obtained explicit approximations (for T = 0). Especially the results 
are valid for the electron gas at low density where plasma oscillations become 
more significant and may have relevance for high-Tc superconductivity. For 
the electron density of conduction electrons in metals the approximation 
obtained gave an energy correction of roughly — leV. 
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For some higher density of electrons comparison is made with an "exact" 
result. It is noted that the Casimir energy is closely associated with the 
RPA for which the statistical mechanical approach will indicate new ways to 
perform modifications to obtain accurate results. 
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